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^ . Abstract. We study baryon asymmetry generation originated from the leptogenesis in the 

Q _ presence of hypermagnetic fields in the early Universe plasma before the electroweak phase 

• transition (EWPT). For the simplest Chern-Simons (CS) wave configuration of hypermag- 

netic field we find the baryon asymmetry growth when the hypermagnetic field value changes 
due to a^-dynamo and the lepton asymmetry rises due to the Abelian anomaly. We solve the 
corresponding integro-diff'erential equations for the lepton asymmetries describing such self- 
consistent dynamics for lepto- and baryogenesis in the two scenarios : (i) when a primordial 
lepton asymmetry sits in right electrons cr; and (ii) when, in addition to cr, a left lepton 
asymmetry for ei and i^eL arises due to chirality flip reactions provided by inverse Higgs 
decays at the temperatures, T < Trl ~ lOTeV. We find that the baryon asymmetry of the 
Universe (BAU) rises very fast through such leptogenesis, especially, in strong hypermagnetic 
fields. Varying (decreasing) the CS wave number parameter /cq < 10"'^ Tew one can recover 
the observable value of BAU, ?7b ~ 10~^, where /cq = 10~''Tew corresponds to the maximum 
value for CS wave number surviving ohmic dissipation of hypermagnetic field. In the scenario 
(ii) one predicts the essential difference of the lepton numbers of right - and left electrons 
at EWPT time, L^r — L^l ~ (/^e_R — AieL)/7Ew = A/x/Tew — 10~^ that can be used as an 
initial condition for chiral asymmetry after EWPT. 
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1 Introduction 

It is well known that, at the high-temperature symmetric phase of the Standard Model (SM) 
all gauge bosons acquire a "magnetic" mass gap ^ g^T, except for the Abelian gauge field 
associated to the weak hypercharge. Such massless hypercharge field in the hot universe 
plasma appears to be a progenitor of the Maxwellian field which evolves after ElectroWeak 
Phase Transition (EWPT). The hypercharge field, in its turn, may arise from phase transi- 
tions in the very early Universe, before EWPT, such as during the infiationary epoch [1]. 

In the absence of hypermagnetic fields the baryon asymmetry of the Universe (BAU) 
can be produced through a leptogenesis. In particular, such a case was considered in Ref. [2] 
assuming that an initial BAU, having a negligible value in our situation, is preserved since it 
is stored in right electrons cr. Then violation of the lepton number due to Abelian anomaly 
in a strong external hypercharge field provides the growth of lepton number from that value. 
Meanwhile lepton- and baryon number evolution proceeds preserving B — L = const. Thereby 
accounting for hypermagnetic fields one can enhance BAU as well. This allows us to assert 
(see also in Ref. [3]) that fermion number "sits" in hypermagnetic field. 

One can easily understand why the authors of Ref. [2] considered the scenario of BAU 
generation with one lepton generation chosen as br. They found that the tiny Yukawa 
coupling to the Higgs for right electrons, he = ^/2me/v = 2.94 x 10~^, provides for such 
leptons the latest entering the equilibrium with left particles through Higgs decays and inverse 
decays. This means that any primordially-generated lepton number that occurred as may 
not be transformed into soon enough switching on sphaleron interactions which wipe out 
the remaining BAU. 

While other right-handed lepton (and all quark) species enter the equilibrium with 
corresponding left particles even earlier, e.g., because of h^^T » he for leptons, they are 
beyond game for BAU generation if their partial asymmetries are zero from the beginning. 

On the other hand, in the presence of a nonzero right electron chemical potential, 
^^eR 7^ 0, the Chern-Simons (CS) term ~ /XeijEyY arises in the effective Lagrangian for the 
hypercharge field [3-5]. It modifies the Maxwell equation in SM with parity violation 
producing additional pseudovector current J5 ~ /ie/jBy in a plasma. It is this current 
which leads to the important ay-helicity parameter in a modified Faraday equation. This 
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parameter appears to be scalar. We remind that the standard magnetohydro dynamic (MHD) 
parameter, which is generated by vortices in plasma, omho ~ — (v-(Vx v))/3, is pseudoscalar. 
Obviously in the isotropic early Universe such vortices are absent, at least, at large scales we 
consider here. 

The Oy'dynamo amplification of a large-scale hypermagnetic field [6] for changing chem- 
ical potential HeR(t) as well as its growth dtfJ-eR > due to the Abelian anomaly in the 
self-consistent hypermagnetic field By were never explored before in literature because of 
the difficulties to solve the corresponding integro-differential equations. In the present work 
we solve that problem in two scenarios. 

First, in the case described in Sec. 2 as saturation regime, we can avoid the adiabatic 
approach BY{t) const and HeR{t) ~ const Ref. [3], using exact solutions of the integro- 
differential equations for dt^eR 7^ (Sec. 3). 

At temperatures T < Trl — 10 TeV chirality flip reactions enter the equilibrium since 
the rate of chirality flip processes, Frl ~ T, becomes faster than the Hubble expansion 
H , Frl > H. This motivates us to consider in the second scenario the extended 

equilibrium state at T < Trl when left leptons enter the equilibrium with cr through inverse 
Higgs decays and acquire non-zero asymmetries ~ HeL{T) = fiu^j^{T) / (Sec. 4). 

The left lepton asymmetries can grow from a negligible (even zero) value at T < Trl 
due to the corresponding Abelian anomaly which has the opposite sign relatively to the sign 
of the anomaly for cr and because left leptons have different coupling constant g'Yi^/2 with 
hypercharge field. Such a difference guarantees the presence of leptogenesis in hypermagnetic 
fields even below Trl all the way down to Tew- Hence it supports generation of the BAU. 

Note that for T > Trlj before left leptons enter the equilibrium with right electrons, 
the anomaly for them was not efficient since the left electron (neutrino) asymmetry was zero, 
fieL = fJ-u^L ~ 0; while a non-zero primordial right electron asymmetry, //gR 7^ 0, kept the 
baryon asymmetry at the necessary level. In other words, for T > Trl the Abelian anomaly 
for left particles was present at the stochastic level, with {Sj^) = = (EyBy) valid only on 
large scales. 

Since fieL = (J-u^l 7^ at temperatures T < Trl there appear additional macroscopic 
pseudovector currents in plasma J5 ~ /ie^Bv which modify ay^helicity parameter governing 
Maxwell equations for hypermagnetic (hyperelectric) fields. This occurs due to the same 
polarization effect described for in Ref. [5] which led to the appearance of CS term in the 
effective SM Lagrangian. 

The kinetics of leptogenesis at temperatures Trl > T > Tew is considered in Sec. 5, 
where we give some details how to solve the corresponding integro-differential equations. 

In Sec. 6 we discuss our results illustrated in Figs. 1 and 2 for the particular case of the 
CS wave configuration of hypermagnetic field. 

2 Saturated regime of baryogenesis 

Let us follow the scenario in Ref. [3] to discern approximations used there and in our approach. 
There are five chemical potentials in symmetric phase corresponding to the five conserved 
numbers. For the three global charges rii = B/3 — Li conserved in SM there are three chemical 
potentials fii, i = 1,2,3, and fiy is the fourth chemical potential which corresponds to the 
conserved hypercharge commuting with the SM Hamiltonian (global {Y) = 0). The fifth 
chemical potential fieR 7^ is the single partial chemical potential for fermions corresponding 
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to the right electron number perturbatively conserved unless Abelian anomaly is taken into 
account. 

Note that formally other partial fermion chemical potentials are zero since there are no 
asymmetries for left leptons including neutrinos, no partial quark asymmetries, etc. There 
are relations between all five chemical potentials (see Eq. (2.16) in Ref. [3]) given by the 
conserved global charges including electroneutrality and hypercharge neutrality of plasma 
(Q) = {Y) = 0. 

Under such equilibrium conditions the last CS anomaly term in the effective Lagrangian 
for the hypercharge field [3], 

L = Vy^^ - J^Y^^ + ^^(By • Y), (2.1) 

originated by the one-loop diagram in FTFT [4] at the finite fermion density {^eR 7^ 0) in hot 
plasma bath of the early Universe (T > Tew) ^ leads to the appearance of parity violation 
terms in the Lagrange equations for hypermagnetic and hyperelectric fields. In particular, 
the modified Maxwell equation takes the form, 

^^"^ + V X By = f J + 4^By1 = cTeond [Ey + V X By + oyBy] , (2.2) 



dt ' V 47r2 

where the last pseudovector term originated by the CS term in Eq. (2.1) is given by the 
hypermagnetic helicity coefficent 

""^^^ = 4.Vcond(T)' ^'-'^ 

that is scalar. In the Lagrangian (2.1) and in Maxwell equation (2.2) the coefficient g' = 
e/ cos 9w is the SM coupling, = (Jo, J) is the vector (ohmic) current with zero time 
component due to the electro-neutrality of the plasma as a whole, Jq = {Q) = 0. In Eq. (2.2) 
we also used the Ohm law J/o"cond = Ey+V x By and introduced conductivity cJcond ~ 100 T 
for the hot universe plasma. 

Note that a more transparent way suggested in Ref. [5] to get the CS term in Eq. (2.1) 
allows to clarify its physical sense as the polarization effect in an external hypermagnetic 
field. 

Using the standard MHD approach one can neglect the displacement current d^y/dt. 
Then we obtain from the Maxwell equation (2.2) the hyperelectric field Ey, 

V X By . . 

Ey = - V X By H - ayBy. (2.4) 

"^cond 

Now one can easily find the change of the CS number density ncs = ^'^(By • Y)/87r2, 



^We changed sign ahead the anomaly term in the SM Lagrangian comparing with what authors had in 
Refs. [5, 6]. This is because of the definition of 75 matrix for right fermions in Ref. [3] which the authors [5, 6] 
rehed on. Such choice of signs in Refs. [3, 5, 6] corresponds to the definition ipu = (1 — 75)i/)/2. Note that 
the change 75 —75 concerns also sign in the Abelian anomaly in Eq. (2.7). After such changes the latter 
(Abelian anomaly) becomes exactly like in the book [8] (pg. 249) where tpR = (1 + j5)tp/2. 



-3- 



that, in its turn, allows us to reveal the generation of BAU, rjB = (?^B — ^b)/'*' hyper- 
magnetic fields, 

3^ A„ _ Z"*""^ ^ By) • By „2 



VBitEw) = - - A?ics = tV / ^ - «y54 dt. (2.6) 

Here s = 27r^r^(7*/45 is the entropy density; g* = 106.75 is the number of relativistic degrees 
of freedom and we substituted in Eq. (2.5) the hyperelectric field from Eq. (2.4). 

So far we neither used the dependence HeRit) on time nor exploited leptogenesis through 
hypercharge fields. To realize this idea authors of Ref. [3] suggested to account for the Abelian 
anomaly for right electrons ^ 

that violates lepton number. In the uniform universe accounting for the right electron 
(positron) asymmetry UeR — ngR = fieRT'^/Q and adding the chirality flip processes (in- 
verse decays crcl — (p^^\ e-R^^eL — ^ ^^~^) with the equivalent rates, Frl and T = 2rpi,L, one 
finds from Eq. (2.7): 

d^^eR 6ff^^(VxBY)-BY 

^ - 4vr2rV,o„d (^B + n^eii, (2.8) 

where we substituted the hyperelectric field from Eq. (2.4) and neglected left electron (neu- 
trino) abundance rising through inverse Higgs decays. Direct Higgs decays do not contribute 
since one assumes for simplicity the zero Higgs boson asymmetry, n^p = n(p. 
The rate of all inverse Higgs decay processes Frl [2], 

r,,.5.3.10-3.^(^)^T.(J^)(l^). (...) 

vanishes just at EWPT time, x = 1. Here he = 2.94 x 10~^ is the Yukawa coupling for 
electrons, Fq = 121, variable x = t/tEW = (Tew/T)'^ is given by the Friedman law, mQ(T) = 
2DT^{1 — Tg^/T^) is the temperature dependent effective Higgs mass at zero momentum 
and zero Higgs vacuum expectation value. The coefficient 2D ~ 0.377 for m'^{T) is given by 
the known masses of gauge bosons ?n^, mw-, the top quark mass rrit and a still problematic 
zero-temperature Higgs mass estimated as mn ~ 125 GeV (see Ref. [9]). 

Note that the helicity term in Eq. (2.4) proportional to aY ~ /^ei? entered Eq. (2.8) as 
— Fb ^eR with the rate 

"^cond 

Then in the adiabatic approximation s = T = and dfieR/dt = 0, we reproduce from 
Eq. (2.8) the right electron chemical potential that is similar to Eq. (2.26) in Ref. [3] (with 
the opposite sign!) 

6g'2(v X By) • By 



^J'eR 



47r2rVeond(F + FB)' 



■^The opposite (here positive) sign in the Abelian anomaly comparing with what was used in Ref. [3] and 
in Refs. [5, 6] corresponds to the standard definition of 75-matrix as in the case of analogous Adler anomaly 
in QED calculated in the book [8]. See our footnote 1 above. 
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Substituting ay with such a chemical potential into the integrand for the baryon asymmetry 
in Eq. (2.6) we reproduce the result, cf. Eq. (2.32) in Ref. [3], 



35' 



,/2 



to 



Tb \ (V X By) • By 



r + FB 



dt 



'^cond 



3g^^r[(VxBY)-BY] 
47r2s(r + rB)acond 



Mo 



(2.10) 



independently of the sign for 75 we used here as it should be. In Eq. (2.10) we put r/B(io) = 
and used a saturated value -By independent of time. Then assuming approximately constant 
integrand we substituted the expansion time ^ew = -Mo/2T^^, Mq = Mpl/1.66^/g*. 

In the scenario of Ref. [3], for the simplest CS wave configuration of the hypercharge 
field, Yo = Y^ = 0, = Y{t) sin k^z, Yy = Y{t) cos koz, for which (V x By) • By = koB^{t), 
where Byit) = koY{t), with the use of the rates T and Fb defined above the integral in first 
line in Eq. (2.10) takes the form, 



rjB{x = 1) =7.3 X 10" 



ko 



IO-^Tew. 

(1 - x')x''^dx' 



( By{x' 



XQ 



[l-x' + 0.16x'2(-By(x')/1020 G)2] V 1020 G 



(2.11) 



Here the parameter A;o/10 ''Tew < 1 is given by the maximum wave number surviving ohmic 
dissipation of hypermagnetic fields at the finite conductivity cJcond — 

100 T. 

For completeness we present here the a2-dynamo formula from Ref. [5] for amplification 



of the hypermagnetic field in the Fourier representation, BY(k,t) = By^ exp[Jl^^^{^aY{t')k 
r}Y{t')k'^ )]dt', with an arbitrary scale A = k~^ = kt/y/cky and k > 1, 



By{x) 



B. 



(0) 



exp 



10"4 



(2.12) 



Here ^e_R = fJ-eR/T is the dimensionless right electron chemical potential, rfy = (dcond)"^ is 
the magnetic diffusion coefficient, and x = t/tEW < 1 is the dimensionless time used above. 
One can see from Eq. (2.11) that varying the two parameters: either the seed hypermagnetic 
field By^ in dynamo formula (2.12) or the CS wave number ko (or both) one can obtain the 
known BAU mitEw) - 6 x 10"^°. 



3 Dynamical evolution of BAU in hypermagnetic fields beyond saturated 
regime 

For the lepton (baryon) asymmetry densities La = Ua — ^^a {B = tlb — n-^) using Hooft's 
conservation law for the electron generation, 

dLeL dLeR _ 1 d(nB - nB)/s _ 1 dr^B , . 

dt dt ~ 3 dt ~ 3 dt ' ^ ■ ^ 

assuming, as in Ref. [3], the absence of left lepton asymmetries, L^l = 0, and the direct 
consequence of the Abelian anomaly for right electrons (2.7), dLgfi/dt = (5''^/47r2s)(EY -By), 
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we get BAU at the EWPT time coinciding with Eq. (2.6) in the adiabatic case s = 0, 

3^/2 ftEw / (V X By) • By „<, \ dt 



EWJ 



47r2 



to 



(VxBv):Bv_^^^, 

<^cond / S 



(3.2) 



Let us use the notation x = t/tsw = (^Ew/^)^ taken from Friedman law and substitute 

LeR = (eRTyes = VrT^ {6s X 10^). 

Again using the CS wave configuration with (V x By) • By = kQBy^x) and relying on 
the total kinetic equation for the right electron asymmetry in Eq. (2.8) where contrary to 
Eq. (2.11) obtained from Eq. (2.8) in the saturation regime dt^eRit) = 0, we get BAU in the 
case dtHeR{t) / 0, 



r]B{x) = 1.14 X 10"^ / dx' 



XQ 



dynix') {l-x') , ■ 

^^ + ro^^yR(x) 



(3.3) 



To obtain Eq. (3.3) we have solved the kinetic equation (2.8) rewritten in new dimensionless 
notations in the form 



dx 



Bqx^/^ - A 



oVR 



B 



(0) 
Y 



where 



Bn = 6.4 



1020 G 



a,3/2gV'{^) 



To ^?/R, 



(3.4) 



^0 = 19.4, 



IO-^Tew, 

are constants chosen for hypermagnetic fields normalized the exponent e^ is given by the 
hypermagnetic field squared, 



By{x) 



B 



and we substituted the hypermagnetic field in Eq. (2.12) with the fastest growth, k = 2, 

yl{x')dx' 



By{x) = B^'^ exp 



10.5 



(3.5) 



that corresponds to the hypermagnetic field scale A = 2r]Y / cxy ■ 

The most important step in the solution of the complicated integro-differential equa- 
tion (3.4) is the following. Accounting for the time derivative de'^^^^ /dx = {21y^^(x) / y/x)e^'^^^ 
with the exponent e'^'^^^ given by Eq. (3.4) itself, or by the first derivative dy-^/dx, let us 
differentiate Eq. (3.4). One obtains the non-linear differential equation of the second order 
which we solve numerically. 



dV_21y| 
dx^ -v/x 



dyR 
dx 



:i-x) 



{dy^/dx + Tq{1- x)yii/Vx) 



x3/2(i?oxi/2_^oyj,) 
X [iBox-lAox'/'yn-Aox'/'^^ 
(1 - x) d^R (1+^,, _n 

+ 1 7= 1 ^ ^ nfn VR — 0. 



^/x dx 



2x3/2 



(3.6) 



^In numerical estimates we substitute the parameter fco/(10 ^ Tew) = 1 that is the upper hmit for the CS 
wave number, ko < 10~^Tew, to avoid ohmic dissipation of hypermagnetic field. 



-6- 





(a) 



(b) 




(c) 




(d) 



Figure 1. Chern-Simons wave configuration for hypcrmagnctic field -By(0 with the maximum wave 
number fco/(10~^ Tew) = 1 surviving ohmic dissipation, (a) Normalized chemical potential t/r = 
10'* X ^eij versus time, (b) Normalized hypermagnetic field versus time, (c) Baryon asymmetry versus 
time. The solid lines correspond to By^ ~ 10^" G, the dashed lines to By "* = 10^* G, and the 
dash-dotted lines to By ■* — 10^^ G. (d) Baryon asymmetry versus time for the small wave number 

(0) 



fco/(10-7 Tew) = 10"^ and = 10'" G 



The initial conditions at xq = (Tew/^rl)^ = 10 ^, Tew 
chosen as yR(xo) = 10~^ (or S,eR{xo) = 10"^'') and 

2 



100 GeV, Trl = lOTeV, are 



dx 



[Bo^/x^ - Aoy^iixo)] 



B 



(0) 
Y 



3/2 



(1 - Xo) 



XQ 



(3.7) 



1020 G / " 

The growth of BAU rj-Q{t) given by Eq. (3.3) is illustrated in Fig. 1(c). 

4 Abelian (triangle) anomaly for left electrons and left neutrinos 

In contrast to QED where the electron lepton number is conserved and Abelian anomaly 
terms in an external electromagnetic field F^y [8], 



dj 



167r2^^" ' 



dj 



dxi^ 



+ 



IGvr 



_7? 7?Mi' 
2 t^'^ ' 
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do not contradict to this law, d^j'^ = 9^0'^ + i^) = 0, the leptogenesis in presence of an 
external hypercharge field Y^i, always exists, i9^(j^ 0; since the coupling constants 

e — )• g'Yji^L/2 are different for right singlet eji and left doublet L = {ueL, ei)'^- 

In addition to Abelian anomaly for right electrons below we assume the presence of 
Abelian anomalies for the left doublet at the stage of chirality flip processes, T < Trl [5], 



^x^' 647r2 



V^^" = -:;^(Ey-By), yL = -i. (4.1) 



These anomalies provide the appearance of the left electron (neutrino) asymmetries. This 
means that the new CS terms in effective Lagrangian with finite /ie^ = I-l^^^ ^ should arise 
due to the polarization mechanism [5]. 

Namely the 3-vector components for the statistically averaged left electron and 
neutrino pseudocurrents. 



Jj5 = [(e7i75e) + {l^eLljlbl'eL)] , 

arise in plasma as additional macroscopic currents modifying Maxwell equations. 



J5 = By. (4.2) 

Here the U(l) gauge coupling fiig') = g'Yh/'^ plays the role of an "electric" charge associated 
to C/y(1), with 1l = —1 being hypercharge of the left-handed electron (neutrino). 

Obviously, /.tgL 7^ HeR in the scenarios of Refs. [2, 3], where /igL = is kept forever, 
while we assume that only as an initial condition, HeL{to) = starting at Tq = Trl- 

If we accept our ansatz with Abelian anomalies for left particles (4.1) we obtain from 
the same Hooft's rule (3.1) the reduced coefficient in the expression similar to Eq. (3.2) 
35'V47r2 3(5'Vvr2)[l/4 - 2 x (1/16)] = 3g''^/8Tr^ and from Eq. (4.2) the modified helicity 
coefficient a™"'^ in the hyperelectric field (2.4), 

Svr^cJcond 

Finally the baryon asymmetry resulting from (3.1) takes the form: 

3g'' f^^w / (VxBy)-By ^modRsV* U 



to \ "^cond / S 



Note that due to Higgs boson inverse decays the chemical potentials Henit) and ^ehii) en- 
tering to oy (q!y°*^) changed over time infiuencing also the amplitude of hypermagnetic 
field Byit) accordingly dynamo Eq. (2.12). 

5 Kinetics of the leptogenesis 

Note that we consider the stage of chirality fiip processes when left electrons and left neutrinos 
supplied by the Higgs decays (inversed decays) in the presence of right electrons enter the 
equilibrium with the same densities at the temperatures below T < Trl , nL = UeL = n-i^eL ' 
having the same chemical potentials /ieL = fJ-u^L- 
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Accounting for the Abelian anomalies for eji and left doublet L = {vcLt^l)'^ we get 
the following system of kinetic equations for the charged lepton asymmetry densities, nji = 

hl/s 



neR — ngR and ul = n^L — ngL, or for the corresponding lepton numbers, LeL 
and LeR = ur/s, [5]: 



L 



dL, 



eR 



r,'2 



dt 
dt 
dt 



^ Ey • By + 2T^i^{LeL - LeR), for e^eL f^^\eRUeL 



47r2 



{-) 



J2 



J2 



9' 



Ey • By + T-Eih{LeR - Lei), for crcl f 
Ey • By + rRL(LeR - Lel), for eRVeL ^ ^ 



(0) 



(+) 



(5.1) 



Here the factor=2 in the first line takes into account the equivalent reaction branches and 
we neglected Higgs boson asymmetries, = n^, or the Higgs decays into leptons do not 
contribute in kinetic equations (5.1). 

One can easily see from the sum of equations (5.1) that in correspondence with BAU 
given by Eq. (4.4) the inverse Higgs processes in the Hooft's rule d?7B/dt = 3[dLe_R/dt + 
dLei/dt + dLj^^^/dt] contribute via indirect way, through the hypermagnetic fields only. 
Thus, baryon number sits in hypermagnetic field itself. 

Below we write differential equations describing evolution of the right and left asym- 
metries that follows from the kinetic equations (5.1). We solved them numerically with the 
corresponding initial conditions and the solutions are illustrated in Fig. 2. 

For right electrons cr we use the notation of the asymmetry yR(x) = 10^ x CeR{x) and 
get from the first equation in (5.1) after differentiation, cf. Eq. (3.6), 



d^VK 21(yR + yL/2)^ 



dx^ 



dx 



(l-x) 



[dy^/dx + (ro(l - x)/^/x){yR - ^l)] 



+ ro 



x3/2[SoxV2_ Ao(yR + 2/L/2 
2^°^ - 2 

(l-x) ( dyR dyL 
dx 



^Bqx - ^AQX^/^{y^ + 2/l/2) - A^x^/^ 



dx 



(1 + x) 

ro^^(yR-yL) = o. 



dyR ^ 1 dyL 
dx 2 dx 



(5.2) 



For left leptons (both and Veh) using the notation yh{x) 
from the second equation in (5.1), 



10'* X £,eL{x) one obtains 



d^yL , 21(yR + yi^/2f 



dx2 



4^ 



dx 



(l-x) 



(yR - yh) 



+ 



1 [dyR/dx + (ro(l - x)/Vi)(yR - yL)] 
4 x3/2[BoxV2_Ao(yR + yL/2)] 

X (^Box - ^A^x^/\y^ + yL/2) - Aqx^/^ 
(l-x) / dyL dyR\ „ fl 



dyR ^ 1 dyL 
dx 2 dx 



2^/x \ dx 



dx 



4x3/2 



(yL - yn) = 0. 



(5.3) 
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The initial conditions at xq = (Tew/^rl)^ = 10 ^, Tew = 100 GeV, Trl = lOTeV, are 
chosen as yR(xo) = 10"^ (or £,eR{xo) = lO"^'') and 



dj/R 



dx 



[Bo^/x^ - AoyB^ixo)] 



B 



(0) 
Y 



3/2 ^ (1 - Xq) 



XQ 



1020 G 



For left particles we choose yh{xo) = S,eL{xo) = and 



(5.4) 



dx 



-T [Bq^/x^ - Aoy^{xo)] 



B 



(0) 
Y 



1020 G 



^3/2 _^ ^°^ j/R(xo). 



(5.5) 



The BAU (4.4) for arbitrary x takes the form, 



, r dyR(x') , ^ (1 - x') 



vb{x) = 0.57 X 10-4 J dx' I + To^^^ [yR(^0 - yL(x')] 



(5.6) 



Accounting for the modified helicity parameter a™"'^ given by (4.3) the a2_(iynamo 
formula for the hypermagnetic field By{x) in the case of the fastest growth, see Eq. (3.5), 
takes the form, 



By{x) = B^Kxp 



10.5 



[?/R(^') + 2/L(x')/2]2dx'^ 



XO 



(5.7) 



6 Discussion 

We study BAU evolution in the electroweak plasma of the early Universe in the epoch before 
EWPT in the presence of a large-scale hypermagnetic field By and assuming a primordial 
right electron asymmetry 7^ that grows in a hypercharge field due to the Abelian 
anomaly. We considered dynamical BAU evolution in the two scenarios: (i) neglecting left 
lepton asymmetry, fi^L = Ul = (Sec. 3); and (ii) assuming its appearance, HeL 7^ 0, at 
temperatures below the chirality fiip processes enter equilibrium, T < Trl ~ lOTeV (Sees. 4 
and 5). 

In the first scenario such a violation of the lepton number given by Eq. (2.7) leads to 
the growth of BAU shown in Fig 1(c) up to the big value ?7b ~ 10"^ — 10"^. We obtained 
that estimate for the largest CS wave number ko = 10"'' Tew surviving ohmic dissipation 
in the case of simplest CS wave configuration for By. Note that decreasing /cq down to 
ko = (IQ-^^ — 10-^4)Tew; we can get the known BAU value r/B — 6 x 10"^'^ for large scales 
of hypermagnetic field, cf. Fig. 1(d). 

The maximum right electron asymmetry at the EWPT time, yR(tEw) — 0.3, is not 
sufficient to amplify the seed hypermagnetic field By^ too much. One can easily estimate from 
the dynamo formula (3.5) that even in the case the fastest growth the hypermagnetic field 
By{x) rises a few times only, cf. Fig 1(a). This is a consequence of our choice of the simplest 
CS wave configuration of By. We think that for the 3D-configuration of hypermagnetic field 
the dynamo amplification should be stronger n2-times, where n = ±1, ±2, ... is the number 
of knots of the tangled By-field. Moreover, for the chosen CS wave configuration a seed field 
By^ should be strong from the beginning to influence the growth of the lepton asymmetry 
T/R and BAU, rye- The stronger By^ the more efficient leptogenesis is (see in Fig. 1). 
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Figure 2. Chern-Simons wave configuration for hypcrmagnctic field -By(0 with tlie maximum wave 
number fco/(10~^ Tew) = 1 surviving ohmic dissipation, (a) Normalized ehemical potential t/r = 
10'* X ^ei? versus time, (b) Normalized chemical potential j/l = 10** x versus time, (c) Normalized 
hypermagnetic field versus time, (d) Baryon asymmetry versus time. The solid lines correspond to 
= 10^0 G, the dashed lines to = lO^i G, and the dash-dotted hues to B^'^ = lO^^ G. 

We conclude that for non-helical CS-wave one needs a preliminary amplification of By^ 
in epoch of inflation to provide both dynamo mechanism and baryogenesis in our causal 
scenario in radiation era. 

In the more realistic second scenario involving left electrons and neutrinos z/gL we 
get the growth of their asymmetries yhix) due to the corresponding Abelian anomalies that 
become efficient at temperatures T < Trl when inverse Higgs decays provide chirality flip 
eR — )• L = {i'eL,eL)'^ (Sec. 4). We reveal the additional CS term in effective Lagrangian 
for the hypercharge field arising due to polarization mechanism [5] that modifies the 
hypermagnetic helicity coefficient in Eq. (4.3) combined from the right and left electron 
asymmetries. As a result the growth of yR{x) and r]-B{x) reduces comparing with the first 
scenario for fieR / alone (see in Fig. 2). 

The important issue in such scenario is a non-zero difference of lepton numbers, L(,r — 
LeL 7^ at the EWPT time, that can be used as a possible starting value for chiral anomaly 
which provides evolution of Maxwellian fields down to temperatures T ~ lOMeV [10]. The 



corresponding difference of asymmetries seen in our Fig. 2(a,b) is estimated at the level 
(,eR — CeL =| -^^/Tew |— 10~^ for = 10^^ G that is a bit below values Ajj/Tew seen in 
Fig. 1 in Ref. [10]. 

We missed here the sphaleron wash-out of BAU due to the appearance of left lepton 
asymmetries believing as in Ref. [2] that time before EWPT during cooling Trl > T > Tew 
is short to evolve such process. We plan to study that effect in future as well as to generalize 
our approach for the case of SD-configuration of hypermagnetic field. 
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There are two kinds of mistakes in our previous work [1]. The values of hypermagnetic 
field in eqs. (3.5) and (5.7) coming from eq. (2.12) were not appropriate for the case of 
Chern-Simons (CS) wave and the values of baryon asymmetry in eqs. (3.3) and (5.6) were 
overestimated due to the incorrect coefficient. 

The a^-dynamo formula (2.12) is applicable for an arbitrary ID-configuration of hyper- 
magnetic fields with the wave number k ~ ay/rjy that, in turn, changes over time through 
the lepton asymmetry, k{t) ~ aylt) ~ Hej^{t). However, assuming such simplest hypercharge 
field configuration as the CS wave, Yq = Yz = 0, = Y{t) s'mkoz, Yy = Y [t) cosk^z, with 
the fixed wave number /cq = const, we should use the exact solution of Faraday equation 



BY{t) = exp l^^o /*[«Y(i') - fcor/Y(i')]di'| , 



:i) 



instead of eq. (2.12). Therefore instead of eq. (3.5) in ref. [1], substituting ay given by 
eq. (2.3), one gets 



BY{t) = b'^^ exp <! 3.5 



IO-^Tew 



0.1 



vr 



kn 



IO-^Tew 



and instead of eq. (5.7) substituting 0^°^ given by eq. (4.3) we obtain 



BY{t) =B^^ exp I 3.5 



kn 



X 



ao-^TEwy 

{yR{x') + yL{x')/2) 
vr 



0.1 



ko 



IO-^Tew 



dx' 



dx' 



(2) 



(3) 



where x = t/tEW < 1 is the dimensionless time and ur^l = 10^(/^eRL/^) = IO^^erl 
lepton asymmetries. One can see from eqs. (2) and (3) that dynamo amplification occurs 
negligible even for the maximum wave number /cq = 10~^rEw for which hypermagnetic field 
survives against ohmic dissipation. Thus, figures 1(b) and 2(c) in ref. [1] for By/B^\ based 
on incorrect eqs. (3.5) and (5.7), were false and we do not replot them here. As a result, 
other plots are almost insensitive to the value ByK Using eq. (2) and differentiating eq. (3.4) 
we get, instead of eq. (3.6), the following non-linear second order differential equation, which 
we solve numerically. 



dx^ 



So (3.48 X lO^^yR - 1.7 x 10"^ x BQ^/^ 



2Box 



0. 



(dyn/dx + ro(l - x)yR/^) 



xV^iBox'/^-Aoyn) 
[1 - x) dyR 



dj/R ^ p (1 
dx 



-VR 



-Aox^/'^yR - Aqx 



3/2 



dyR 
dx 



X 



dx 



r (1 + ^) 



2x3/2 



(4) 



The initial conditions for this equation are correct including eq. (3.7). 
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Using eq. (3) we get, instead of eqs. (5.2) and (5.3), 

^2 



-^-Bo(3A8 X 10-1 X (yR + yL/2) - 1.7 x 10"^ x Bq^ 

[dyn/dx + (ro(l - x)/^/x){yR - y^)] 
x3/2[SoxV2_Ao(yR + yL/2)] 

2^02; - ^Aox^/\yn + yL/2) - Aqx^/^ 



d^R _^ 1 dyL 
dx 2 dx 



0" 



(1 



dj/R _ dj/L \ 
dx dx J 



(1 + x) , 

ro^;^(yR-yL) 



0, 



(5) 



dV^^Bo 
dx2 4 

X 



3.48 X 10-1 X (yR + yi^/2) - 1.7 x 10' 
x) 



X Bo^/x 



d2/R , p (1 
+ 1 0— 



+ 



dx <y^-y^\ 

1 [dyK/dx + (ro(l - x)/V^)(j/R - t/l)] 

4 x3/2[i?oxl/2_ylo(yj, + yL/2)] 

f 2Sox - ^^o^'/'(yR + yL/2) - ^0X^/2 



dj/R ^ IdyL 
dx 2 dx 



(1 - x) ( dyL dyR 
2^^ 



dx dx 



^ (1 + x) , , 
ro^^(yL-yR) 



0, 



(6) 



obeying the same initial conditions (5.4) and (5.5). We show in figure 2(a,b) the evolution 
of lepton asymmetries yR,L(a^) given by solutions of eqs. (5) and (6) instead of figures 2(a,b) 
in ref. [1]. 

We correct also wrong values of the baryon asymmetries (BAU) in eqs. (3.3) and (5.6) in 
ref. [1], where the number of relativistic degrees of freedom, g* = 106.75, was not accounted 
for in the denominator. It led to the incorrect plots in figure l(c,d) and in figure 2(d) 
correspondingly. The correct expressions read 



?7b(x) = 1.068 X 10' 



dx' 



f dyR(xO 
\ dx' 



(1 



X 



(7) 



instead of eq. (3.3) in scenario with right electron asymmetry yR = 10"^ x (/Xgj^/T) evolution 
and 



7?b(x) = 5.34 X 10" 



' I f dyR(xO {I-X') r . , 



(8) 



dx' vx 

instead of eq. (5.6) in scenario when both lepton asymmetries yR and yL = 10^ x (/igj^/T) 
evolve. 

Below using the numerical solution of eq. (4) for yR(x) shown here in figure 1(a) we 
replot figure 1(c) in ref. [1] as figure 1(b) given by eq. (7) for the same maximum wave number 
of the Chern-Simons wave = 10"'^r. We replot also figure 1(d) in ref. [1] as new figure 1(c) 
given by eq. (7) changing = 10"i^r used in ref. [1] to k^ = IQ-^^T. Such wave number 
provides observable BAU at t = tEW) ~ 10" i*^. 

We also present BAU evolution in figure 2(c,d) instead of figure 2(d) in ref. [1] using 
here the correct eq. (8) based on the solution of eqs. (5) and (6) for asymmetries yR,L(2;). The 
curves in figure 2(c) are plotted for the same /cq = 10"^T and in figure 2(d) for k^ = lO-^^T. 
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Figure 1. (a) Normalized chemical potential ?/r = 10^ x versus time, (b) Baryon asymmetry 
versus time. Panels (a) and (b) are built for fco/(10~^ Tew) = 1- The solid lines correspond to 
B'y^ ^ 10^" G, the dashed lines to = lO^i G, and the dash-dotted hues to S^"' ^ lO^^ G. (d) 
Baryon asymmetry versus time for the small wave number ko/{10~'^ Tew) = 10"^ and B^^ = lO^o G. 
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Figure 2. (a) Normalized chemical potential j/r = 10^ x versus time, (b) Normalized chemical 
potential ?/l ~ 10"^ x ^p,^ versus time, (c) Baryon asymmetry versus time. Panels (a)-(c) arc built for 
ko/{10~'^ Tew) = 1- The solid hues correspond to B^'' = lO^o G, the dashed lines to B^'' = lO^i G, 
and the dash-dotted lines to By ^ = 10^^ G. (d) Baryon asymmetry versus time for the small wave 
number ka/ilQ-'^ Tew) = 10"^ and S^"^ = lO^o G. 
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